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Abstract 

Charge Coupled Devices (CCDs) have been successfully used in several high energy physics 
experiments over the past two decades. Their high spatial resolution and thin sensitive layers 
make them an excellent tool for studying short-lived particles. The Linear Collider Flavour 
Identification (LCFI) collaboration is developing Column-Parallel CCDs (CPCCDs) for the 
vertex detector of a future Linear Collider. The CPCCDs can be read out many times faster 
than standard CCDs, significantly increasing their operating speed. A test stand for measuring 
the charge transfer inefficiency (CTI) of a prototype CPCCD has been set up. Studies of the 
CTI have been performed at a range of readout frequencies and operating temperatures. 
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Abstract — Charge Coupled Devices (CCDs) have been suc- 
cessfully used in several high energy physics experiments over 
the past two decades. Their high spatial resolution and thin 
sensitive layers make them an excellent tool for studying short- 
lived particles. The Linear Collider Flavour Identification (LCFI) 
collaboration is developing Column-Parallel CCDs (CPCCDs) for 
the vertex detector of a future Linear Collider. The CPCCDs can 
be read out many times faster than standard CCDs, significantly 
increasing their operating speed. A test stand for measuring the 
charge transfer inefficiency (CTI) of a prototype CPCCD has 
been set up. Studies of the CTI have been performed at a range 
of readout frequencies and operating temperatures. 



L Introduction 

The study of radiation hardness is crucial for the application 
of CCD detectors in High Energy Physics experiments [l]-[3]. 
The LCFI collaboration has been developing and testing new 
CCD detectors for about 10 years [l]-[5]. Previous experimen- 
tal results on CCD radiation hardness were reported for exam- 
ple in [6]-[8]. Several models increased the understanding of 
radiation damage effects in CCDs [9]-[ll]. The measurements 
and analyses reported in this paper have been carried out 
in the LCFI collaboration [3]. Simulation and modelling of 
CCD radiation hardness effects for a CCD prototype with 
sequential readout was reported at IEEE '2005 [12], comparing 
full TCAD simulations with analytic models was reported at 
IEEE' 2006 [13] and in Ref. [14], simulation and modelling of 
a CCD prototype with column parallel readout (CPCCD) at 
IEEE'2007 [15] and in Ref. [16]. 

This work focuses on experimental measurements and a 
method to determine the charge transfer inefficiency (CTI) 
performed with the CPCCD at a test stand at Liverpool 
University. 

The environment with high radiation near the interaction 
point at a future Linear Collider creates damage to the CCD 
material which leads to defects acting as electron traps in the 
silicon, as indicated in Fig.[T] The mechanism of creating traps 
has been discussed in the literature, for example in Refs. [17]- 
[19]. These traps result in charge transfer inefficiency. In a 
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Fig. 1. Schematic view of how radiation damage is created in phosphorus- 
doped silicon. Two types of vacancy defect are created: Vacancy-Oxygen 
(electron trap at 0.17 eV below the conduction band) and Vacancy-Phosphorus 
(electron trap at 0.44 eV below the conduction band). 

phosphorus-doped device, two types of traps are created. The 
first one is relatively shallow with energy level 0.17 eV below 
the conduction band and the second is deep with energy level 
0.44 eV below the conduction band. These traps result in 
charge transfer inefficiency. 

The column parallel technology is in development to cope 
with the required readout rate. CPCl is a two-phase CCD 
prototype capable of 50 MHz readout frequency. In this paper 
we demonstrate the method to determine the CTI value with 
an un-irradiated CPCCD (CPCl). 

II. Test stand for CCD operation and cryostat 

UNIT 

A test stand has been set up with readout electronics 
and a cryostat unit. The temperature range of the cryostat 
is from room temperature down to about —140 °C. This 
temperature has been achieved with cold nitrogen gas by 
boiling liquid nitrogen. The very low operating temperature 
is required to obtain sensitivity to the CTI peak structure for 
0.17 eV traps where the peak position is expected to be near 
— 130 °C [16]. Initial measurements have been performed on 
an un-irradiated device in standalone mode, where the signal 
from four columns of the CCD were amplified and connected 
to external ADCs. An ^^Fe source emitting 6 keV X-rays was 
attached to a holder at a distance of 5 cm from the CCD to 
provide the signal charge. Figure [2] shows a picture of CPCl 




Fig. 2. CPCl (in the center of the picture) with external electronics. 

and the associated control and front-end readout electronics, 
which were also placed in the cryostat, wherefrom the four 
outputs were fed to rack-mounted amplifiers and ADCs. 

It was observed that the CTI increased sharply when the 
amplitude of the sinusoidal clock pulses applied to read out 
the CCD was reduced to less than 2 V peak-to-peak. In order 
to suppress any effect from the clock amplitude, settings for 
each data acquisition were tuned to produce 3.00 Vpp clock 
pulses. The performance of the drivers limited the maximum 
clocking frequency to about 22 MHz. The settings depended 
on the frequency and had to be adjusted for each cryostat 
temperature. 

III. Signal measurement and removing charge 

SHARING BY USING A 3x3 CLUSTER METHOD 

The fast ADCs convert the signal charge after amplification 
with a wideband preamplifier. The four columns are read out 
in 4 channels by ADCs. Three columns (channels 1 to 3) are 
adjacent whereas channel 4 is separate. The signal charges of 
400 pixels plus 10 overclocks were acquired in 3000-5000 
frames per measurement. This leads to an accurate statistical 
precision. The collected data has been analysed using the 
ROOT package [20]. First, we begin by applying correlated 
double sampling, where the difference between the signals 
of two consecutive pixels is taken to be the signal charge 
collected by the latter pixel. This decreases electronic noise. 
Next, we remove common mode noise components for each 
column with the help of the overclocks. Overclocks are charges 
read out after the last pixel of the column. We average all 
overclocks of a column and subtract the average from the data 
of all pixels in the column. As an example for channel 2 Fig. |3] 
shows the distribution of ADC codes for all recorded pixels 
and all frames of one measurement. The analysis is performed 
on column 2 because it is situated between columns 1 and 3. 
This allows to apply a 3 x 3 cluster method. 

Figure [3] illustrates that there are 3 regions in ADC codes. 
The first region (I) has a high number of pixels with low 
ADC codes from electronics noise only (no charge from X- 
rays). The noise is fitted by a Gaussian function in order to 
determine the noise threshold. The charge- sharing region (II) 
represents the signal charge shared between 2 or more pixels. 
Region (III) is the signal region where the charge is collected 
in a single pixel that will be analysed to determine the CTI. 
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Fig. 3. Distribution of ADC codes for channel (column) two. Three regions 
are observed: (I) the high peak region which represents the noise, (II) the 
region separating the two peaks which represents the charge sharing between 
pixels, and (III) the X-ray peak region which represents the collected charge 
in a single pixel. The noise peak is fitted by a Gaussian function to determine 
the noise threshold. The data was taken at —80 °C with 8 MHz readout 
frequency. 
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Fig. 4. Illustration of the 3 x 3 cluster method. A signal charge (X-ray) is 
only accepted if no X-ray charge is present in the neighboring pixels. The 
X-ray charge is defined by its threshold determined after a first fit of the X-ray 
signal peak with a Gaussian function. 



In order to remove hits with shared charges, a 3x3 cluster 
method is used (Fig. [4]). This method is only applied to the 3 
adjacent columns (1,2 and 3). Column 2 is considered, where 
the signal charge of a pixel is accepted only if no X-ray charge 
is present in the neighboring pixels. 

The noise threshold is used to separate noise from shared 
X-ray charges. It is determined by fitting the noise peak with 
a Gaussian function. The threshold is usually set to be So- 
above the fitted peak centroid. After removal of the noise and 
reduction of the charge- sharing events, a tail remains for the 
X-ray peak as shown in Fig. |5] In order to separate the X-ray 
signal region from the charge- sharing region we exclude all 
ADC codes below an X-ray threshold which is determined by 
fitting the X-ray peak. We usually set the threshold to be 2a 
below the fitted peak centroid but with weaker X-ray peaks 
the fitted Gaussian sometimes failed to follow the real shape 
of the peak, and then the limits were set by visual inspection. 
In this way we obtain the measured X-ray signal data in the 
pixels with full X-ray charge (Fig. |6]). 
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Fig. 5. X-ray signal peak with a tail from some charge sharing even after 
applying the 3 x 3 cluster method. The distribution is fitted with a Gaussian 
function to determine the threshold to be used to completely remove the charge 
sharing. This is a different dataset than that in Fig. [3] which was obtained 
using a higher preamplifier gain. 



Fig. 7. Distribution of average ADC code in the noise region versus the pixel 
number at 2 MHz and — 60 °C with an exponential fit. The noise region was 
delimited by a 3cr noise threshold. Fit: amplitude A = 0.117 ± 0.056 and 
exponent B = (-5.29 ± 2.85) x 10-^. 
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Fig. 6. X-ray signal peak with no charge sharing. Fig- 8. Distribution of average ADC code in the noise region versus the 

pixel number at 4 MHz and — 6 °C with an exponential fit. The noise region 
was delimited by a 3a noise threshold. Fit: amplitude A = 15.64 ±0.44 and 
exponent B = (-16.24 ± 1.33) x 10"^. 



IV. Baseline correction for different operating 

TEMPERATURES AND READOUT FREQUENCIES 

Owing to an RC coupling at the preamplifier input, the 
output voltage decays exponentially with a time constant equal 
to RC (trc is of the order of 100 /isQ This is expected 
to give rise to an exponential baseline that varies during the 
clocking sequence, adding a contribution to the readout signal 
that depends on the pixel number. The shape of the baseline 
can be determined by comparing the average output level in 
the pixels to each other. For each pixel an average of the 
ADC codes was computed excluding the frames were the 
code exceeded the noise threshold. The average ADC code 
can usually be expected to be an exponential function of the 
pixel number with the RC decay. Figures [7] [Sj [9] and [10 



^For a readout frequency of 4 MHz reading out the whole column of 400 
pixels takes a time comparable to trc- 



show the distributions of the average ADC codes for all 
frames of a measurement as a function of the pixel number 
for different readout frequencies and operating temperatures. 
These distributions are fitted by an exponential function given 
by Aex.p{—Bj), where A corresponds to the signal charge at 
the first pixel, B is the slope and j is the pixel number. This 
function will be subtracted from the X-ray data. Figures [7] 
show approximately zero baseline level, whereas 
and |9] show non-zero baseline fits. The baseline 
level is much smaller at low temperatures, for example at 
-60 °C (Fig. [7]) and at -109 °C (Fig. [To]). Figures [8] and [9] 
show baselines for higher temperatures — 6 °C and —17 °C, 
respectively. The observed baselines show different frequency 
and temperature dependence from what is expected from the 
RC effect. This indicates that other electronic effects also 
influence the baseline. 
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Fig. 9. Distribution of average ADC code in the noise region versus the pixel 
number at 8 MHz and — 17 °C with an exponential fit. The noise region was 
delimited by a 3cr noise threshold. Fit: amplitude A = 8.31 ±0.31 and 
exponent B = (-11.86 ± 1.74) x 10"^. 
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Fig. 11. Fit to average ADC codes without baseline removal at 2 MHz and 
-31 °C. CTI = (2.99 ± 2.69) x 10"^. 




Fig. 10. Distribution of average ADC code in the noise region versus the Fig. 12. Fit to events (hits) without baseline removal at 2 MHz and -31 °C. 

pixel number at 8 MHz and —109 °C with an exponential fit. The noise CTI = (3. 69 ±1.31) x 10~^. 
region was delimited by a 3a noise threshold. Fit: amplitude A = (0.082 ± 
5.035) X 10-3 and exponent B = (1.073 ± 16.853) x 10"^. 



V. Determination of Charge Transfer Inefficiency -g 3«)p 



The charge transfer inefficiency (CTI) in one pixel is defined 
as the ratio of signal lost during transfer (captured by traps) to 
the initial signal charge. For an un-irradiated CCD we expect 
that the CTI value is consistent with zero within uncertainties. 
In order to determine the CTI we first make the overclock 
correction and apply the 3 x 3 cluster method. Then, the X-ray 
peak was fitted with the Gaussian function to determine the X- 
ray threshold. This was used to construct the two distributions 



with a fit of average ADC codes as shown in Fig. 11 and with 
a fit of individual events (X-ray hits) as shown in Fig. [12] Both 
distributions are fitted with the first-order polynomial function 
Po + Pij, where Pq corresponds to the charge at the first pixel, 
Pi is the slope and j is the pixel number. These distributions 



are plotted without baseline removal. Figures [13] and [14] show 
the distributions of the averages and all events as a function of 
pixel number with baseline removal. The CTI is determined 
for the four cases using CTI = —Pi/Pq. 
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Fig. 13. Fit to average ADC codes with baseline removal at 2 MHz and 
-31 °C. CTI = (-0.09 ± 2.70) x 10-^. 
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Fig. 14. Fit to events (hits) with baseline removal at 2 MHz and 
CTI = (0.63 ±1.33) X 10-5. 
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Fig. 15. CTI versus temperature at 2 MHz readout frequency. The CTI is 
shown resulting from fits to average ADC codes and to events (hits), with 
and without baseline removal. The shown error bars result from the precision 
of the fits. 
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Fig. 16. CTI versus temperature at 4 MHz readout frequency. The CTI is 
shown resulting from fits to average ADC codes and to events (hits) , with 
and without baseline removal. The shown error bars result from the precision 
of the fits. 
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VI. CTI RESULTS PRE-IRRADIATION FOR DIFFERENT 
READOUT FREQUENCIES 



Figures 15 16 and 17 show the CTI values as a function 



of temperature for an un-irradiated CPCl at different readout 
frequencies. The CTI has been calculated using two methods, 
a linear fit of averages versus pixel number and using a linear 
fit of events (X-ray hits) versus pixel number. As it is expected 
for an un-irradiated CCD, the CTI is small because the density 
of traps is small. In the temperature region below —25 °C no 
frequency dependence of the CTI measurements is observed. 
One may note that in the region of rather high temperatures 
(above —25 °C) the effect of baseline removal is large and 
the determined CTI values indicate that the baseline effect is 
overestimated. The effect of the baseline removal decreases 
with increasing frequency as the readout clocking sequence 



becomes shorter compared to trc (as visible in Fig. 17 for 
8 MHz readout frequency). 



Fig. 17. CTI versus temperature at 8 MHz readout frequency. The CTI is 
shown resulting from fits to average ADC codes and to events (hits) , with 
and without baseline removal. The shown error bars result from the precision 
of the fits. 



VII. Conclusions and Outlook 

An un-irradiated CPCCD is operated in a range of temper- 
atures from —10 °C to —136 °C (liquid nitrogen cooling) 
with different readout frequencies 2, 4 and 8 MHz. The 
spectrum of a ^^Fe source is measured with this device. 
The CTI is analysed for different readout frequencies and 
operating temperatures. A clear X-ray signal is extracted by 
identifying isolated hits (3x3 method). The baseline is 
subtracted. The CTI value is small and compatible with zero 
as can be expected for an un-irradiated CPCCD. Further CTI 
measurements with a CPCCD after irradiation and refinement 
of the analysis method are planned. 
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